THE ASTROPHYSICAL JOURNAL LETTERS, 941:L28 (4pp), 2022 December 20 
© 2022. The Author(s). Published by the American Astronomical Society. 


OPEN ACCESS 


https: / /doi.org / 10.3847 /2041-8213 /aca8a0 


CrossMark 


Interstellar Meteors Are Outliers in Material Strength 


Amir Siraj® and Abraham Loeb® 
Department of Astronomy, Harvard University, 60 Garden Street, Cambridge, MA 02138, USA; amir.siraj@cfa.harvard.edu, aloeb@cfa.harvard.edu 


Received 2022 September 20; revised 2022 December 1; accepted 2022 December 4; published 2022 December 16 


Abstract 


The first interstellar meteor larger than dust was detected by US government sensors in 2014, identified as an 
interstellar object candidate in 2019, and confirmed by the Department of Defense in 2022. Here, we describe an 
additional interstellar object candidate in the CNEOS fireball catalog and compare the implied material strength of 
the two objects, referred to here as IM1 and IM2, respectively. IM1 and IM2 are ranked first and third in terms of 
material strength out of all 273 fireballs in the CNEOS catalog. Fitting a log-normal distribution to material 
strengths of objects in the CNEOS catalog, IM1 and IM2 are outliers at the levels of 3.50 and 2.60, respectively. 
The random sampling and Gaussian probabilities, respectively, of picking two objects with such high material 
strength from the CNEOS catalog are ~1074 and ~1076. If IM2 is confirmed, this implies that interstellar meteors 
come from a population with material strength characteristically higher than meteors originating from within the 
solar system. Additionally, we find that if the two objects are representative of a background population on random 
trajectories, their combined detections imply that ~40% of all refractory elements are locked in meter-scale 
interstellar objects. Such a high abundance seemingly defies a planetary system origin. 


Unified Astronomy Thesaurus concepts: Interstellar objects (52) 


1. Introduction 


CNEOS' on 2014 January 8, a meteor detected by the US 
Department of Defense (DOD) sensors through the light that it 
emitted as it burned up in the Earth’s atmosphere off of the 
coast of Papua New Guinea in 2014, was determined to be an 
interstellar object in 2019 (Siraj & Loeb 2019), a conclusion 
that was confirmed by independent analysis conducted by the 
DOD in 2022 (Shaw 2022) although the numerical uncertain- 
ties were not provided to the scientific community. The object, 
which we refer to here as IM1, predated the interstellar object 
‘Oumuamua by 3.8 yr and the interstellar object Borisov by 5.6 
yr. The measured peak flare apparent in the light curve of IM1 
at an altitude of 18.7 km implies ambient ram pressure of ~194 
MPa when the meteor disintegrated (Siraj & Loeb 2022a). This 
level of material strength is >20 times higher than stony 
meteorites and >2 times larger than iron meteorites. IM1 was 
also dynamically unusual—its speed relative to the local 
standard of rest (LSR) is shared by less than 5% of all stars. 

In this Letter, we describe an interstellar meteor candidate 
from the CNEOS catalog mentioned in Table 8 of Pefia- 
Asensio et al. (2022) and in Section 5 of Siraj & Loeb (2022b), 
which we refer to as IM2. We then explore the statistical 
likelihood that interstellar meteors reflect the same distribution 
of material strength as noninterstellar meteors. 


2. Interstellar Meteor Candidate 


The Python code implemented here used the open-source N- 
body integrator software REBOUND? to trace the motion of the 
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meteor under the gravitational influence of the solar system 
(Rein & Liu 2012). 

We initialize the simulation with the Sun, the eight planets, 
and the meteor, with geocentric velocity vector (VXobs, VYobs» 
VZobs) = (—15.3, 25.8, —20.8) kms ', located at 40.5° N 
18.0° W, at an altitude of 23.0 km, at the time of impact 
t; = 2017 March 9, 04:16:37 UTC, as reported in the CNEOS 
catalog. We then use the IAS15 adaptive time-step integrator to 
trace the meteor’s motion back in time (Rein & Spiegel 2015). 
This does not account for air drag, which would lead to an even 
higher impact speed and therefore heliocentric speed, given the 
encounter geometry. The slowdown of IM1 due to air drag was 
estimated in earlier work (Siraj & Loeb 2022a). 

There are no substantial gravitational interactions between the 
meteor and any planet other than Earth for any trajectory within the 
reported errors. Based on the geocentric impact speed reported by 
CNEOS, vobs = 36.5kms', and the heliocentric impact speed 
was ~50kms'. We find that the meteor was unbound with an 
asymptotic speed of v, ~ 25.9 kms_' outside of the solar system. 

We find the heliocentric orbital elements of the meteor at 


impact to be the following: semimajor axis a=-—1.1 au; 
eccentricity e = 1.6; inclination i = 26°; longitude of the 
ascending node Q = —12°; argument of periapsis w = 241°; and 


true anomaly f = 300°. IM2 was vsr ~ 40 kms! away from 
the velocity of the local standard of rest (Sch6nrich et al. 2010). 

Given the explosion energy of ~4 x 10" erg and the 
atmospheric impact speed of ~36.5kms~', we adopt equiva- 
lence between the preexplosion kinetic energy and the energy 
in the explosion, finding that the objects mass was 
~6.3 x 10°g. A comparison between the properties of IM1 
and IM2 is included in Table 1. 

We note that the fact that both IM1 and IM2 have low orbital 
inclinations (i < 30°) is puzzling since interstellar objects are 
expected to have a uniform distribution in cos i. Specifically, the 
random likelihood of two orbital inclinations of i < 30° drawn 
from a uniform distribution in cosi is (1 — cos(30°))? ~ 2%. 
However, 2I/Borisov had an inclination of 44°, and the 
likelihood of drawing three inclinations of i<45° is 
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Table 1 
Comparison between IM1 and IM2 
Designation ti Vobs View VLSR Y; m n 
(km s7!) (km s~!) (km s7’) (MPa) (g) (au) 
IM1 2014-1-08 17:05:34 44.8 42.1 60 194 4.6 x 10° 18 x 10° 
IM2 2017-3-09 04:16:37 36.5 25.9 40 75 6.3 x 10° 2.7 x 10° 
(1 — cos(45°))? ~ 3%, so there may be an inclination bias in Fit Values: 0.65 and 0.47 
the source of a certain class of interstellar objects. 11/ 1.24 a 
‘Oumuamua had an orbital inclination of i=122°, more 
indicative of a background population uniform distribution. 2 TA 
Q 
g 
= 
3. Material Strength Comparison a cal 
As a meteor travels through the atmosphere, it experiences z 0.6. 
friction due to air. The dynamical pressure is pv*, and the 5 
dynamical pressure corresponding to the peak power in the J 
meteor light curve describes the material strength of the meteor g Sii 
since crossing a certain level of ram pressure causes the object = 
to deform and break apart (Trigo-Rodriguez & Llorca 2006, 0.2 
2007; Popova et al. 2011). 
Based on estimates for comets and carbonaceous, stony, and 0.0 


iron meteorites (Chyba et al. 1993; Scotti & Melosh 1993; 
Svetsov et al. 1995; Petrovic 2001), Collins et al. (2005) 
established an empirical strength—density relation for impactor 
density p; in the range 1-8 g cm *. The upper end of this range 
gives a yield strength of Y;~ 50 MPa, corresponding to the 
strongest known class of meteorites, iron (Petrovic 2001). Iron 
meteorites are rare in the solar system, making up only ~5% of 
modern falls (Zolensky et al. 2006). 

We computed the ram pressure at breakup for all 273 
fireballs in the CNEOS catalog. Interestingly, IM1 and IM2 
display the first- and third-highest material strengths, respec- 
tively, among all of the fireballs. Figure | is a histogram 
showing all of the fireballs in the catalog and highlighting IM1 
and IM2. Figure 2 shows the light curves for IM1 and IM2 with 
the peak ram pressures highlighted. 

The probability of randomly drawing two of the top three 
material strengths out of all 273 fireballs is ~(3/273) ~ 10+. 
Therefore, if IM2 is confirmed to be an interstellar meteor, simple 
random drawing dictates that there would be a ~99.99% chance 
that interstellar meteors come from a population with material 
strength characteristically higher than meteors originating from 
within the solar system, a notion suggested by Pefia-Asensio et al. 
(2022). 

The material strength data follow a log-normal distribution, 
shown in Figure 1 with mean p = 0.65 and o = 0.47 MPa. As a 
result, IM1 and IM2 represent 3.50 and 2.60 deviations from 
the mean, respectively. These deviations correspond to 
2.4 x 1074 and 4.5 x 107° single-tailed probabilities, respec- 
tively. Combining these independent events, we find a ~10~° 
probability of getting the material strengths of IM1 and IM2 by 
random chance. This Gaussian perspective implies a 
~99.9999% chance that interstellar meteors are characteristi- 
cally stronger than meteors from within the solar system. 

We note that the ram pressure could be overestimated by the 
fact that DOD satellites only detected the brightest sections of their 
luminous paths. For instance, the US government sensor data on 
the Chelyabinsk event gives a dynamic strength 2 or 3 times as 
large as that measured in recovered meteorites (Pefia-Asensio et al. 
2022). 


-1.0 -05 00 05 10 15 20 25 
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Figure 1. Histogram of ram pressure at breakup pv’ for all fireballs in the 
CNEOS catalog with altitudes and velocities reported at the times of peak 
brightness. Navy blue and bright blue correspond to IMI and IM2, 
respectively, which are ranked first and third in terms of material strength 
out of all 273 fireballs. The black line indicates the best fit to a Gaussian with 
mean and standard deviation matching the data 4. = 0.47 and ø = 0.65. 
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Figure 2. Total power released in the IM1 and IM2 fireballs (navy blue and 
bright blue lines, respectively) as a function of ram pressure pv’. Peak 
brightness is reached at 194 and 75 MPa, respectively, for the two fireballs. 
Typical stony and iron meteorite yield strengths 1-5 and 50 MPa, respectively, 
are indicated for convenience of comparison. Note that 1 TW = 10” ergs |, 
and 1 MPa = 10’ dynes cm~?. The IM2 light curve is calibrated based on the 
total energy released (taking account of a normalization error in the published 
light curve). There is an early flare in the IM2 light curve indicating some 
amount of lower-strength material in addition to the clearly central flare 
corresponding to a metallic composition. 


4. Implications for Local Mass Budget 


For a background population on random trajectories drawn 
from an isotropic distribution in the LSR, the number density 
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Figure 3. IM1 and IM2 in size-abundance parameter space, expressed as 
number per star per differential unit of log size. The vertical error bars 
correspond to 95% Poisson uncertainties, while the horizontal error bars 
correspond to a factor of 2 in mass in each direction. 


implied by the detection of an interstellar meteor is 


T 
+ 1 
i VooTR ELL + (Vesc/Voo)7] 9 


where I is the implied rate, va is the speed outside of the solar 
system,” Rọ is the radius of the Earth, and vesc = ./2GM>/ de 
is the escape speed from the Earth’s position, where Mo is the 
mass of the Sun and d, the distance between the Earth and the 
Sun. Gravitational focusing is accounted for by the term in the 
denominator [1 + (Vese/Voo)7]. We adopt T ~ 0.1 yr for both 
IM1 and IM2 (Siraj & Loeb 2019) and speeds outside of the 
solar system of vo ~ 42kms™' and v,,~26kms_', respec- 
tively. We find that the number density implied in the 
detections of IM1 and IM2 are nmı ~ 1.8 x 10° au~? and 
nyo ~ 2.7 x 10° au~3, as shown in Figure 3. 

Given the respective masses of ~4.6 x 10° g and ~6.3 x 10° g, 
we find that the detections of IM1 and IM2 imply, respectively, 
ambient local abundances of ~1.2 Mẹ pe * and ~25 Mẹ pe > of 
similar objects. 

The local stellar mass density is ~0.04 Mo pe ° (Bovy 2017). 
The local density of the interstellar medium is 1.2 cm? (McKee 
et al. 2015), implying ~0.03 Mo pe. All refractory elements 
(metals and silicates) sum to a total mass fraction of ~0.3% at 
solar metallicity, implying that the local budget of metals and 
silicates in stars and dust is ~70 Mẹ pe’. We conservatively 
assume that IM1 and IM2 are composed of refractory elements 
even though their material strengths imply that they were 
primarily metallic in composition. If IM2 is indeed an interstellar 
object, the detections of IM1 and IM2 together imply that ~40% 
of all refractory elements locked from stars and the ISM are 
locked in meter-scale interstellar objects. Refractory elements are 
also locked in objects bound to the Sun, with the largest 
theorized reservoir being the Oort cloud (Stern & Shull 1990). 


3 Without knowing the velocity dispersion of the population, our best 
estimate is taking the velocity of the object relative to the Sun as a “typical 
relative speed,” including the velocity of the Sun relative to the LSR and the 
velocity dispersion of the population in the LSR. We only have a sample of two 
objects, which are not sufficient for a better statistical analysis to find the 
velocity dispersion of the population. 
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5. Discussion 


If interstellar meteors are formed in planetary systems, the 
natural limit to the scale of mass ejected is the total budget of 
the minimum mass solar nebula model (MMSN), which is of 
order ~1% of stellar mass (Weidenschilling 1977; Hayashi 
1981; Desch 2007; Crida 2009). The result reached here 
indicates that if IM2 is confirmed an interstellar object, the 
detections of IM1 and IM2 combined imply that ~2/3 of the 
mass budget in stars is necessary to provide the refractory 
elements to produce a population of interstellar meteors that 
would make the detections of IM1 and IM2 likely. This result 
thereby provides a new constraint on planetary system 
formation since it requires nearly 2 orders of magnitude more 
mass than the MMSN (Weidenschilling 1977; Hayashi 1981; 
Desch 2007; Crida 2009). Note that the mass budget discussed 
exceeds that of objects “‘Oumuamua-sized and larger,” which 
itself is an unsolved puzzle (Loeb 2022; Siraj & Loeb 2021). 
The extraordinary mass budget required to produce interstellar 
meteors seemingly defies planetary system origins and suggests 
some other highly efficient route for creating meter-scale 
objects made of refractory elements. Interestingly, there is a 
paucity of refractory elements observed in the gas phase in the 
interstellar medium (Savage & Sembach 1996; Maas et al. 
2005; Delgado Inglada et al. 2009), an observation that could 
potentially reflect refractory elements being locked in inter- 
stellar objects. Supernovae have been observed to produce 
iron-rich “bullets,” which could be a possible origin of IM1 and 
IM2 (Loeb et al. 1994; Strom et al. 1995; Stone et al. 1995; 
Wang & Chevalier 2002; Tsunemi & Katsuda 2006; Perret & 
Timmes 2009; Miceli et al. 2013; Tsebrenko & Soker 2015; 
Sandoval et al. 2021). 
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